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Cosmological models with cold dark matter composed of 
weakly interacting particles predict overly dense cores in the 
centers of galaxies and clusters and an overly large number 
of halos within the Local Group compared to actual obser- 
vations. We propose that the conflict can be resolved if the 
cold dark matter particles are self-interacting with a large 
scattering cross-section but negligible annihilation or dissipa- 
tion. In this scenario, astronomical observations may enable 
us to study dark matter properties that are inaccessible in the 
laboratory. 



Flat cosmological models with a mixture of ordinary 
baryonic matter, cold matter, and cosmological constant 
(or quintessence) and a nearly scale- invariant, adiabatic 
spectrum of density fluctuations are consistent with stan- 
dard inflationary cosmology and provide an excellent fit 
to current observations on large scales 1 Mpc) la How- 
ever, an array of observations on galactic and subgalactic 
scales (< few Mpc) appears to conflict with the struc- 
ture formation predicted by analytical calculations and 
numerical simulations. The predictions are based on the 
standard view of cold dark matter as consisting of parti- 
cles with weak self-interactions, as well as weak interac- 
tions with ordinary matter. 

A generic prediction for weakly self-interacting dark 
matter, independent of other details of the cosmological 
model, is that cold dark matter forms triaxial halos with 
dense cores and significant dense substructures within the 
halo. Yet, lensing observations of clustersB reveal central 
regions (roughly galactic scale) with nearly spherical low 
density cores. Dwarf irregular galaxies appear to have 
low density coresmj with much shallower profiles than 
predicted in numerical simulationsE'cl The persistence 
of bars in high surface brightness galaxies imply that 
galaxies like our own Milky Way also have low density 
cores.el Observations atthe Local Group reveal lessJ-haa. 
one hundred galaxies ,E2L pihile numerical simulations^' t3 
and analytical theorjOEJ predict that there should be 
roughly one thousand discrete dark matter halos within 
the Local Group. 

In this paper, we propose that the inconsistencies with 
the standard picture may be alleviated if the cold dark 
matter is self-interacting with a large scattering cross- 
section but negligible annihilation or dissipation. The 
key feature is that the mean free path should be in the 
range 1 kpc to 1 Mpc at the solar radius, where the dark 
matter density is about 0.4 GeV/cm 3 . The large scatter- 



ing cross-section may be due to strong, short-range in- 
teractions, similar to neutron-neutron scattering at low- 
energies, or weak interactions mediated by the exchange 
of light particles (although not so light as to produce a 
long-range force). Depending on the interaction and the 
mean free path, the requisite mass for the dark matter is 
in the range 1 MeV to 10 GeV. For the purposes of our 
proposal, only two-body scattering effects are important 
so either repulsive or attractive interactions are possible. 
Exchanged particles should be massive enough that they 
are not radiated by the scattering of dark matter parti- 
cles in the halo. 

We are led to consider self-interactions because ordi- 
nary astrophysical processes are unlikely to resolve the 
problems with standard, weakly interacting dark matter. 
Consider the dwarf galaxy-problem. One might suppose 
that supernova explosionst^l could cause the galactic core 
density to be made smoother; but, while the explosions 
suppress-star formation in dwarf galaxies, numerical sim- 
ulationalj find that starbursts in dwarfs are very ineffi- 
cient at removing gas or matter from the core. One might 
also consider whether the apparent overabundance of ha- 
los found in simulations can be explained if the low ve- 
locity halos form primarily low surface brightness galax- 
ies ,E3 which are difficult to find. However, while low 
brightness galaxy surveys suggest a steeper luminosity 
function outside of groups J13 even these surveys do not 
find enough small galaxies to eliminate the discrepancy 
between theory and observations. If star formation in 
dwarfs is sufficiently suppressedjlj then they should have 
been detected as gas clouds in the local groupcEl or ex- 
ternal systems. HI surveys-.do not find large numbers of 
small isolated gas clouds. tiJ Even if any of the processes 
were successful in reducing the number of visible dwarfs, 
the dense small halos would still persist. When these ha- 
los fall onto galactic, djftks, they will heat the stellar disks 
and destroy themOE3'E2l These dense halos will also set- 
tle to the centers of the central halo and produce a high 
density core in galaxies and clusters. Since the baryon 
fractioxL-in the centers of low surface brightness galaxies 
is lowjLZl hydro dynamic processes are not likely to alter 
their dark matter profiles n u 

The success of the cold dark matter model on large 
scales suggests that a modification of the dark matter 
properties may be the best approach for resolving the 
problems on small scales. If the dark matter is not cold, 
but warm (moderatefv-jjelativistic), this alleviates some 
of these discrepancies. Eil However, the remarkably good 
agreement between standard cold dark matter (CDM) 
models and the observed power spectrum of Lyman a ab- 
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sorbersa likely rules out warm dark matter candidates. 

We propose that a better resolution is dark matter 
that is cold, non-dissipative, but self-interacting. There 
are stringent constraints on the_.mteractions between 
dark matter and ordinary mattercScland on long range 
forces between dark matter particles. E3 However, as long 
as the dark matter annihilation cross-section is much 
smaller than the scattering cross-section, there are rel- 
atively few constraints on short-range dackpLnatter self- 
interactions. Carlson, Machacek & HalSEil suggested 
a self-interacting dark matter model in which the dark 
matter particle is warm rather than cold. Their model 
assumed that the dark matter plus ordinary matter sum 
to the critical density predicted by inflationary cosmol- 
ogy. Their purpose was to reduce the power on 10 Mpc 
scales in the dark matter mass spectrum as required if 
the normalization of the spectrum is to agree fluctua- 
tions measured by the COBE satellite. Subsequently, de 
Laix et alL3 pointed out that the alteration cannot simul- 
taneously fit the IRAS power spectrum and the observed 
properties of galaxies. Our proposal does not suffer from 
this problem because we assume a cosmological model 
with a low matter density and a cosmological constant (or 
quintessence) which satisfies the COBE constraint with- 
out self-interactions. Our proposed self-interactions do 
not change structure on the 10 Mpc scale but only on 
the 1 kpc scale. Consequently, our model satisfies the 
constraints raised by de Laix et al. 

To be more specific, we suggest that the dark matter 
particles should have a mean free path between ~ 1 kpc 
to 1 Mpc at the solar radius in a typical galaxy (mean 
density 0.4 GeV/cm 3 ), for reasons to be explained be- 
low. For a particle of mass m x , this implies an elastic 
scattering cross-section of 
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intriguingly similar to that of an ordinary hadron. (In 
this paper, we consider the case of dark matter particles 
scattering only from themselves but, in a forthcoming 
paper, we consider the possibility that dark matter is 
a stable, neutral hadron.) If the dark matter particles 
scatter through strong interactions similar to low-energy 
neutron-neutron scattering, then the cross-section is a — 
Atto 2 , where a is the scattering length. For neutrons, the 
scattering length is more than 100 times its Compton 
wavelength. Using the estimate a ~ IOO/to^T 1 , we obtain 
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Alternatively, the self-interaction may be weak but 
longer-range, as in the case of the exchange of a light 
intermediate vector boson of mass m y , in which case the 
cross-section is a w a y m 2 /my. The mass of the vector 



boson must be large enough that there is no dissipation 
when dark matter particles scatter; this requires that 
m y > 450 eV (m x / GeV)(v/ 200km /s) 2 , where v is the 
typical velocity of dark matter particles in the halo. This 
mass scale for m y corresponds to a force that is short- 
range compared to the dark matter interparticle spacing 
(about 1 cm in the halo). Hence, we need only consider 
two-body interactions in our analysis. If m y = gm x and 



0(1), then the maximum dark matter mass is 
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Beyond what is expressed in the relations above, there 
is no significant constraint on how light the dark matter 
particles can be. 

The strong self-interaction might occur if the dark mat- 
ter consists of particles with a conserved global charge 
(such as a hidden baryon number) interacting through 
a hidden gauge group (e.g., hidden color). If the gauge 
group is unbroken, then the particles experience strong 
interactions which can be non-dissipative but the par- 
ticle number is conserved. M-theory and superstrings, 
for example, suggest the possibility that dark matter 
fields reside on domain walls with gauge fields sepa- 
ratedJmm ordinary matter by an extra (small) dimen- 
sion.E3E3 Similar scenarios can be constructed in purely 
four-dimensional supergravity models. Note that, if the 
sum of the hidden baryon number and the ordinary sec- 
tor baryon number is zero, then Q x — (m x /m pro ton)^b — 
0.19 (m x /4GeV) (using VL b h 2 = 0.02 and h = 0.65). The 
particles we suggestinclude light versions of WIMPZIL- 
LASB and Q-ballsH 

How does the mean free path of the dark matter par- 
ticles affect astrophysics? Since interactions only alter 
the evolution of cold dark matter when the density inho- 
mogeneities are large, the cosmic microwave background 
(CMB) and large-scale power spectrum measurements 
are not sensitive to the self-interactions. So long as the 
dark matter is cold (T x /m x < $, where $ is the depth of 
the gravitational potential), the dark matter will collapse 
to form a bound halo regardless of its collisional proper- 
ties. If the dark matter mean path were much longer than 
~ 1 Mpc, the typical dark matter particle would not ex- 
perience any interactions as it moves through a halo. In 
this regime, the usual, triaxial cold dark matter halo with 
dense core forms through gravitational collapse. On the 
other hand, if the dark matter mean free path is much 
smaller than 1 kpc, then the dark matter behaves as a 
collisional gas and this alters the halo evolution signifi- 
cantly. The dark matter will shock: this will heat up the 
low entropy material that would usually collapse to form 
a core and produce a shallower density profile. Since col- 
lisions tend to make the dark matter velocity distribution 
isotropic, the halo can not be triaxial and will only be el- 
liptical if flattened by significant rotation. Since dark 
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halos form with little angular momentum, if the dark 
matter is not dissipative, then all halos .will be nearly 
spherical. X-ray observations of clusterscll reveal that 
most halos are moderately ellipsoidal. This implies that 
the collision time scale for dark matter near the half-mass 
radius of clusters must be longer than the Hubble time: 
one of the strongest constraints-cn this model. Studies 
suggest that polar ring galaxiesES are only mildly triax- 
ial and oblate with the equatorial plane of the dark halo 
nearly coinciding with that of the stellar body. If the 
dark matter has an isotropic distribution function and 
the baryons form a disk, then the dark matter will form 
a slightly flattened halo. 

In our scenario, we consider a mean free path in the 
intermediate regime, larger than 1 kpc, but smaller than 
~ 1 Mpc. Particles in this range have 1 — 10 3 interac- 
tions per Hubble time in the local neighborhood, which 
is overdense by 10 6 relative to the mean density of the 
universe. At the virial radius of a typical galactic halo, 
which is overdense by ~ 200 relative to the mean density 
of the universe, the typical particle has less than 1 colli- 
sion per Hubble time. Thus, near the virial radius, halos 
can have anisotropic velocity ellipsoids and will be triax- 
ial. However, in the inner halo of galaxies, dark matter 
is collisional. Infalling dark matter is scattered before 
reaching the center of the galaxy so that the orbit dis- 
tribution is isotropic rather than radial. These collisions 
increase the entropy of the dark matter phase space dis- 
tribution and lead to a dark matter halo profile with a 
shallower density profile. The characteristic scale for the 
core would correspond to an "optical depth" of 1, the 
"photosphere" of the dark matter. 

When a dwarf halo with a low velocity dispersion falls 
into a larger, high velocity dispersion halo, the high ve- 
locity particles will scatter off the low velocity particles. 
After the collision, neither particle is likely to be bound to 
the dwarf. As dark matter is slowly removed, the dwarf 
halo expands, this also makes the halo more vulnerable 
to tidal stripping and shock heating. This process will 
slowly evaporate all substructure in the larger halo, par- 
ticularly, in the centers of galaxies, groups and clusters. 
However, near the half-mass radius, the collision time is 
longer than the Hubble time, so substructure will only be 
destroyed in the inner portions of halos. This dark mat- 
ter evaporation will protect galaxy disks from dynamical 
heating by collisions with dwarfs in the halo. Dwarfs with 
high central densities evaporate more slowly as particles 
at optical depth greater than ~ 1 are shielded by other 
particles from collisions. Intriguingly, all of our Galaxv^s. 
dwarf companions have very high phase space densities. t£l 
As the central density and particle velocities increase, not 
all collisions will lead to particles being deflected out of 
the dwarf halo. A small fraction of the particles will be 
rescattered within the spheroidal core of the dwarf and 
a fraction of the momentum absorbed. (The phase space 
for capture is suppressed by (er/u) 3 , where a is the ve- 



locity dispersion in the dwarf core and v is the typical 
particle velocity in the halo.) This could produce a ram 
pressure drag that can slow the dwarf halo and cause it 
to spiral into the cores of larger halos. Numerical stud- 
ies of this process are required to determine whether the 
observed fraction of denser dwarfs are likely to survive 
after a Hubble time. 

The halos of large (e.g. L*) galaxies moving through 
groups and clusters are less prone to destruction. When a 
cluster dark matter particle strikes a galaxy dark matter 
particle, the probability that the recoiling particles will 
escape from the galaxy is significantly less than unity. 
For cross-sections near the smaller end of our suggested 
range, most collisions take place within the half mass ra- 
dius, where the escape velocity from a galaxy halo is com- 
parable to the characteristic recoil velocity. Thus, many 
collisions will not lead to recoil energies large enough to 
escape. At the smaller end of our cross-section range, 
the probability of a typical particle experiencing a col- 
lision during a Hubble time approaches unity only for 
galaxies that fall deep into the cluster core. Towards the 
larger end of our suggested range, the galaxies are opaque 
to dark matter and the collision products are also likely 
to experience multiple collisions within a massive galaxy 
halo. Because of these multiple collisions, the collision 
products are unlikely to escape from the galaxy. This 
effect is more important in massive halos as the optical 
depth through a galaxy scales as the one-third power of 
its mass, Thus, for large galaxies, the background cluster 
or group will primarily heat the entire dark halo rather 
than evaporate dark matter from it. 

The presence of collisions will lead to energy transport 
within the dark matter halo, which eventually leads to 
core collapse.Ea We can obtain an estimate of the core col- 
lapse time from Quinlan'scJ Fokker-Planck simulations of 
the evolution of an isolated cluster of interacting particles 
with a central density profile of r~ l and an outer profile 
of r~ 3 . These models have a temperature inversion in 
the core and undergo two stages of core collapse. During 
the first stage, the inner region expands as heat is trans- 
ported inwards. After 0.1 half-mass relaxation times, the 
inner 1% of the mass has moved out in radius by factor 2. 
After roughly 3 half- mass relaxation times, the whole sys- 
tem collapses as heat is transported from the virial radius 
outwards. We suspect that this second stage is delayed 
in our cosmological context as the infall of new mate- 
rial is constantly adding heat near the dark halo photo- 
sphere. If we model our Galaxy as starting with a den- 
sity profile corresponding to the Navarro- Frenk- White fit 
to CDM modelsQ with y 2 oo = 225 km s -1 , VL = 0.3, 
Hq = 65 km s _1 Mpc -1 and dimensionless concentration 
parameter c = 8, then its half mass radius is 220 kpc. 
Since the density at the half mass radius is 300 times 
smaller than the local dark matter density, the particle 
mean free path at that radius is between 0.3 - 300 Mpc 
and at that radius, the particle is in the weakly interact- 
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ing regime. Thus, for our Galaxy, the core collapse time 
(roughly 3 half-mass relaxation times) is between 4.5 - 
6000 Gyr, or perhaps significantly longer if the collapse 
stage is delayed by the infall of new material. Hence, for 
most of our range of parameters, the collapse time for our 
Galaxy exceeds the lifetime of the universe, yet there is 
sufficient number of interactions to lower the dark matter 
density in the inner 5 kpc of our Galaxy. As the particle 
mean free path approaches the lower bound (0.3 Mpc) or 
the upper bound (300 Mpc) at the half-mass radius, our 
estimates suggest that one or the other condition is not 
satisfied, but more accurate methods are needed to deter- 
minc the precise ra nge. [Note added in proof: A. Burkert 
(jastro-ph/000240S| ) has constructed an N-body code that 
simulates self-interactions and obtains results consistent 
with these estimates based on the Fokker-Planck approx- 
imation.] 

To summarize, our estimated range of a/m for the 
dark matter is between 0.45-450 cm 2 /g or, equivalcntly, 
8 x 10~( 25 ~ 22 ) cm 2 /GeV. Numerical calculations are es- 
sential for checking our approximations and refining our 
estimates. Even without numerical simulations, we can 
already make a number of predictions for the properties 
of galaxies in a self-interacting dark matter cosmology: 
(1) the centers of halos are spherical; (2) dark matter ha- 
los will have cores; and (3) there are few dwarf galaxies in 
groups but dwarfs persist in lower density environments; 
and, (4) the halos of dwarf galaxies and galaxy halos 
in clusters will have radii smaller than the gravitational 
tidal radius (due to collisional stripping). Intriguingly, 
current observations appear to be consistent with all of 
these predictions. 
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